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ABSTRACT 


Atmospheric turbulence introduces random phase distortions in optical 
imaging systems. The development of new laser and imaging systems requires 
information on the spatial and temporal distribution of this atmospheric 
turbulence. Measurements of the image spread and the jitter induced by the 
atmosphere on an optical system provide two techniques to quantify these 
phenomena. This thesis evaluates a Spectra Sources Lynxx PC Plus charge 
coupled device (CCD) array as an atmospheric turbulence sensor. Data 
acquisition and processing programs were written to measure the image spread 
of a point source and centroid jitter of a point Source imaged through the 
atmosphere. Since atmospheric jitter measurements require high image frame 
rates on the order of 200 images per second, a large portion of this thesis 
involved measurements of the times for the CCD detector, interface board and 
IBM compatible computer to perform their tasks. Recommendations for higher 


performance are presented. 
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| INTRODUCTION 


This thesis considers a method used to measure atmospheric turbulence 
using a silicon charge coupled device (CCD) camera. Light from a distant star 
consists of almost perfect plane waves before entering the earth’s atmosphere. 
When the waves traverse the atmosphere, local variations in the index of 
refraction induce warping of the wave fronts. Therefore, a telescope pointed at 
a star will produce a distorted image, whose spatial and temporal characteristics 
depend upon the properties of the atmospheric path followed by the starlight. 
Typically, the atmosphere introduces both image spread and temporal jitter, so 
that a CCD array placed at the telescope’s focal plane will reveal a highly 
unstable and complex irradiance pattern. It has been found that both a high 
frame rate (> 200 per second) and correspondingly short exposure time (< 5 
ms) are necessary to avoid under-sampling the dynamic variations of typical 
images. 

The goal of this thesis was to find a way to measure atmospheric turbulence 
with a CCD camera by optimizing the speed and efficiency of the camera 
system’s software and hardware. The software was written in the C language. 
Software components that can affect system speed include both the various 
algorithms used and the overall efficiency of the compiled program code. Three 


compilers’ outputs were compared for executable code speed and utility: 


Microsoft Quick C 2.5, Borland C++ 2.0 and Borland Turbo C++ 1.0. Fixed 
hardware components included a Texas Instruments TC211 CCD image sensor 
in a Spectra Source PC Lynxx Plus camera attached to an 8 bit IBM PC 
compatible interface card. The camera system was installed in two IBM 
compatible personal computers, one using an Intel 80386 and the other a 80486 
processor, and the overall speeds were compared. Finally, the camera-computer 
interface and the CCD image sensor were analyzed to determine whether such 
a low-cost commercial system would be suitable for turbulence measurements, 


or whether a customized device would have to be fabricated. 


ll BACKGROUND 


A. ATMOSPHERIC TURBULENCE 

Turbulence in the atmosphere causes many problems. In astronomy, it 
introduces distortion that obscures image detail. Since exoatmospheric use is 
not practical for most telescopes, Knowing the spatial and temporal distribution 
of turbulence is essential in the planning stages of a new facility or in diagnostic 
evaluation of test systems. Figure 1 sketches a wave front incident on a turbulent 
region [Ref. 1]. Random phase fluctuations in the index of refraction field 
produce scintillations and image blurring. Twinkling or scintillation arises from 
the interference of starlight that traverses multiple paths through the atmosphere. 


Image blurring arises from the reduction in spatial coherence of the wave front. 


B. PARAMETERS FOR THE MEASUREMENT OF TURBULENCE 

The atmospheric index of refraction depends on both temperature and 
pressure. Since pressure fluctuations disperse rapidly, at the speed of sound, 
temperature fluctuations are the main contributor to atmospheric optical 
turbulence. Structure functions provide a way to quantify the statistics of 
atmospheric turbulence. For propagation through the atmosphere, optical 


parameters will involve an integral of the refractive structure function along the 
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Figure 1. Optical atmospheric turbulence distorts a coherent phase front. 

optical path. The primary optical parameter is the spatial coherence length, r.. 
lt represents the transverse autocorrelation length of the electromagnetic field. 
A related parameter is the isoplanatic angle, @,. It represents the angle between 
two different paths where the Strehl ratio of an ideal adaptive optics system is 
within e' of perfect correction [Ref. 2]. Another parameter, the Greenwood 
frequency, f,, represents the electrical bandwidth needed to remove atmospheric 


phase distortions with an adaptive optical system. 


1. Structure Function 
Atmospheric turbulence produces localized variations in the index of 
refraction along an optical path. Tatarski [Ref. 3] defines the structure function 


by 
J, UZ aig) = \ bla) ae Gin) 


where < > represents an ensemble spatial average, r, and r, represent the 
location of two points in space, and x represents an atmospheric parameter of 
interest, such as temperature or index of refraction. Assuming the turbulence is 
isotropic, homogeneous and incompressible [Ref. 4] the Kolmogorov theory of 


turbulence [Ref. 5] shows that 


Dea epree oe (2) 


x 
where C.° is a structure parameter and r,, is the distance between the two points 
in space. By combining equations (1) and (2), the thermal structure parameter, 
C,°, which characterizes the mean squared temperature difference between two 


points in space is [Ref. 3], 


where T.-T, is the temperature difference and r,, is the distance between the two 
points in space. There is a similar expression for the index of refraction structure 
parameter, C,*. Taking the partial derivatives of the atmospheric index of 


refraction [Ref. 3], 


! (4) 


where P is the atmospheric pressure in mbar and T is the Kelvin temperature, 
and assuming isobaric turbulence allows us to write C,* in terms of C,* [Ref. 6], 


=o 


2 _ 
Ch = 72 


(5) 


The pressure is assumed to be constant, since small random pressure 


differences disperse rapidly. 


2. Spatial Coherence Length 


The spatial coherence length, r,, as stated above, measures the 


y 
transverse autocorrelation length of a wave front. It represents the effective 
aperture diameter of a diffraction-limited optical system with a similar angular 
resolution as the system under study. Typical values of r, range from two or 
three centimeters for high turbulence to as much as thirty centimeters for low 
turbulence, when measured from the earth’s surface upward. When observing 
a star from the ground with an optical device, the C,* component of r, is 
cumulative along the optical path, even though most of the degradation of the 


wave front occurs close to the ground. Integrating the optical turbulence C,°(z) 


along the path [Ref. 6], 





O 
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where k is the wave number(2r/Aa), @ is the zenith angle, L is the vertical path 


length and W(z) is a weighting function. Typically W(z) has one of three forms: 


for a plane wave, W=1; for a spherical diverging wave, W=(z/L)°*; and for a 


*® The factor of sec @ compensates for the slant 


converging wave, W=(1-2/L) 
path through a horizontally stratified atmosphere. 

Temperature fluctuations along a vertical path have been measured by 
two means: by microthermal probes carried by a meteorological balloon and by 
an echo sounder [Ref. 7 and 8]. Both of these methods provide credible values 
of r,. However, there is a more direct method for calculating the coherence 
length at the earth’s surface using an optical point source irradiance distribution 
[Ref. 9]. The discrete irradiance values from a two dimensional irradiance 
distribution represent the point spread function, P(x,y). The line spread function 
for x, L(x), is the summation of all of the y irradiance values from P(x,y) for each 
x. L(y) is a similar summation of the x values for each y. Then take the Fourier 
transform of L(x) or L(y) to get the optical transfer function (OTF). Assuming the 
turbulence is isotropic, that is symmetric with respect to rotation about the path, 
and that it is laterally stationary, the absolute value of the OTF equals the 
modulation transfer function (MTF). For an optical system imaging a star the 


observed quantity, MTF, is a product of the MTF’s of the source or star, the 


atmosphere and the instrument [Ref. 9], 
MTF, = MTF,‘MTF,‘MTF,. OW) 


MTF, for a star corresponds to a point source at infinity and is unity. MTF, isa 


measured instrumental transfer function. To determine the atmospheric 
component, MITF,, from the observed quantity, MTF,, we divide MTF, by the 
measured instrumental transfer function, MTF.. 

Fried shows that atmospheric MIF, is a function of the atmospheric 
wave structure function, which is an integral of the index of refraction structure 


parameter along the path [Ref 2]. MTF, reduces to a form [Ref. 9], 


-3.44 (ARV) 5/3 
MIF, =e 0 (8) 


where R is the effective focal length, A is the wavelength of light and v is the 
Spatial frequency. To determine the atmospheric parameter r,, we took a one 
dimensional fast Fourier transform of a line spread function of a star image. 
Dividing this by the known instrument function gives the spatial spectrum similar 
to Equation (8). The next step was to determine the e”' frequency for this 


experimental spatial spectrum from which [Ref. 9], 
I = 25 Ae (9) 


where v, was the e” frequency of the MTF,,. 

Yet another way to calculate r, uses image centroid jitter, also called 
beam wander, by calculating the x and y centroid standard deviation from a 
sequence of point source images. It tS equal to aF, where a is the atmospheric 
angle of arrival fluctuations and F is the optical system's effective focal length. 


The angle of arrival is sketched in Figure 2. Mathematically, angle of arrival 
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Figure 2. The angle of arrival is the least-squares plane that fits the electric field 
vector. 


fluctuations are equivalent to least-squares planes that fit the electric field across 
an aperture. In terms of the index of refraction structure parameter [Ref. 5 and 


10}, 


1.026 
p13 





0? = ee; (2) dz, (10) 
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where D is the diameter of the aperture and L is the optical path. Combining 


equation (10) with equation (6) for a plane wave, 
eles Gee ceeie.’ als’ >, (11) 


To calculate r, from the centroid motion requires a very high image frame rate to 
"freeze" the atmosphere. Characteristically at least 200 frames per second are 
required in order to avoid aliasing that would otherwise underestimate the true 


centroid variance. 


3. !soplanatic Angle 
The tsoplanatic angle measures the angular coherence in the vicinity 
of an object. It defines a cone that constrains the angles over which an adaptive 


optical system will provide valid correction. In terms of C,* [Ref. 1], 


ih & (2.91. k? f °C (z) PAIN 0 fedilear aap (12) 
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where the parameters are the same as those of r,. Simple instruments exist to 
measure @, [Ref. 11]. For this reason, only the coherence length is addressed 


in this thesis. 


4. Greenwood Frequency 
The Greenwood frequency was introduced tn an earlier paper by Darry! 
Greenwood [Ref. 12 and 13]. It is a mean temporal frequency of an adaptive 
optics system and depends on the wind velocity, V, and C,* along the path. For 


a Kolmogorov turbulence spectrum, the Greenwood frequency is [Ref. 13], 


jae Dee A if Cane) Ve za oe (13) 
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where A is the wavelength of light, z is the distance along the propagation path, 
V(z) is the wind speed and L is the distance from the source to the receiving 


optical device. 


C. CHARGE COUPLED DEVICE (CCD) 

An atmospheric coherence length sensor needs an imaging detector to 
measure the point spread function. The speed, reliability and availability of 
charge coupled devices (CCD's) makes them worthy of consideration. The basis 
of a CCD is a metal-oxide-semiconductor (MOS) capacitor forming potential wells 
and channels that comprise light-sensitive pixels and read-out registers. Figure 
3 is asketch of the MOS capacitor structure [Ref. 14]. Electrons excited by the 


photoconductive effect are trapped in potential wells formed under the positively 





Metal (transparent) layer (gate) 


S10, 


p-type silicon 


Figure 3. The metal-oxide-semiconductor (MOS) capacitor structure,the basis of 
the charge coupled device (CCD). 


charged metal gate contact. The number of electrons trapped is proportional to 


the integrated irradiance during an exposure. 
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There are several methods to read-out the device. Voltages with different 
phases placed on the gates transfer charge from pixel to pixel or potential well 
to potential well. A set of electrodes connected together is called a phase. 
Figure 4 shows an example of a three phase device with three pairs of gates (G) 


and three lines (L) connecting the electrodes [Ref. 14]. A positive electrical pulse 


hd . 
A eG Ce A Oe 
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Figure 4. Three pairs of gates (G) are connected with the 
lines (L) to form three different phases. 


from a clock forces a transfer of charge from one phase of pixels to the next. 
Clocking a phase of cells reduces the barrier between pixel columns allowing 
charge transfer. A bucket of water running down several steps with boards or 
barriers on each step is a visual representation of this process. The water flows 
down to a lower potential each time the board is lifted up. The number of 
phases used to transfer out charge packets represents the number of steps. 
There are two methods used to insure that efficient charge transfer occurs 
in one direction. The first uses multiple phases to separate the charge packets. 


The Texas Instruments TC211 [Ref. 15] uses the second method which is to force 


12 


an asymmetry in each well using an ion implantation zone between each pixel. 
The clocking method is the monophase mode, also referred to as a 1+ '% phase 
mode, sketched in Figure 5 [Ref. 16]. One phase “ee intermediate voltage 
level while the voltages applied to the other phase vary on both sides of this 
level. 

The charge packets transfer down channels. There are two. types of 
channels, surface and buried. The water analogy also works here. In a buried 
channel, there is less loss of charges, similar to a pipe. A surface channel, 
however, has a higher capacity, like a canal, but this CCD has higher noise from 


surface imperfections. 
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Figure 5. The monophase clocking method holds a phase at an intermediate 
level while the voltages are applied on either side. 
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During the read-out of a CCD, image smear will occur since light will still 
produce charge in the moving CCD rows, unless a shutter is incorporated. 
Different techniques exist to reduce image smear during read-out. Frame transfer 
devices overcome this disadvantage by effectively having two arrays, one for 
image exposure and a second for storage. Moving charges from the image area 
to an opaque, masked set of pixels quickly allows another exposure to 
commence while reading the storage area. Figure 6 shows two examples [Ref. 


14]. 
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Figure 6. Two schemes for CCD read-out. 


frame transfer method. 


Top, the interline method; bottom, 
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ll! DISCUSSION 


A. EQUIPMENT 

This thesis investigated a Spectra Source Lynxx PC Plus CCD Imaging 
System [Ref. 17 and 18]. The original intent was to use this CCD camera with 
a telescope and personal computer to develop software to measure the 
atmospheric coherence length, r,. During this process it became clear that the 
off-the-shelf CCD system could not provide sufficient frame rates to achieve the 
intended goal. Consequently, the task shifted to determine the causes of the 
sluggish acquisition rates. The CCD detector array, computer interface card and 


the IBM compatible PC were investigated. 


1. CCD Camera 
The camera contains a 192 X 165 pixel Texas Instruments TC211 CCD 
image sensor. Its specifications are in Appendix A [Ref. 19]. Figure 7 is a 
photograph of the sensor centered inside the camera. The other components 
inside the camera head translate the clock signals from the interface board and 
amplify the CCD output. 
Figure 8 shows a functional block diagram of a TC211 CCD image 
sensor [Ref. 19]. Key components include the silicon matrix of pixels and a set 
of gates to shift the charge. The image area gate (IAG) performs a parallel shift 


of all rows for each clock pulse. As a row enters the output serial register the 
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Figure 7. Photograph of inside the CCD camera head of the Spectra Sources 
Lynxx PC+ system. 


serial register gate (SRG) shifts the row to the charge detection amplifier. The 
TC211 includes an antiblooming gate (ABG) that should not be used since it 
introduces a severe nonlinear photo response. 

The method used to acquire and to read-out an exposure determines 


the CCD sensor’s read-out rate. The electrons collect in potential wells, or pixels 
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Figure 8. Lynxx PC+ image sensor CCD, Texas Instruments TC211 block 
diagram. 


after exposing the image area to light. The charges then shift down to the output 
register as rows, through 165 timing cycles. The output register then shifts the 
charges 210 times for each row. This cycle of 210 shifts includes 12 dark pixels 
for a dark reference and 6 dummy pixels used to transfer the charges out of the 
register. Figure 9 shows the charge transfer process [Ref. 19]. Another factor 
to consider for an image sensor chip is the noise. The noise equivalent signal 
for the TC211 is 150 electrons. It depends on (kTC)'? /q, where C is the 
capacitance of the read-out charge collector, k is the Boltzmann constant, T is 
the absolute temperature and q is the charge of an electron. 

The integration time needed to achieve a particular signal-to-noise ratio 


depends on both the telescope optics and the object's radiance. It will be a 
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Figure 9. TC211 image area transfer process, uses a virtual phase. 


constant time added to the read-out time regardless of the number of pixels used 
or transfer times. The image area gate (IAG) and serial register gate (SRG) 


maximum pulse rates determine the sensor read-out rate. 
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The maximum rate for the IAG parallel row transfer gate is 1.5 MHz. 
This corresponds to 0.667 microseconds (us) per row. Clocking the IAG shifts 
the whole array down by one row. This can be used to clear the CCD. To clear 
all 165 rows requires only 110 ps. Since an IAG clock pulse clears the serial 
register, the IAG should not be clocked any faster than the SRG can read-out a 
row of pixels. 

The maximum rate for the SRG single row shift gate is 10 MHz. This 
corresponds to 0.1 us per pixel. Figure 8 shows the pixel read-out order. Six 
dummy pixels will appear before any valid data and the twelve trailing dark pixels 
can be ignored. At the maximum clock rates, the ideal minimum time required 
to read an entire row of pixels is 210 pixels X 0.1 ps = 21 ps. The minimum time 
needed for all 165 rows is 165 rows X 21 ps = 3465 us. 

The total read-out time must include the serial shift time and two 
parallel array shifts. One set of IAG parallel shifts are needed to move each row 
into the output shift register and a second set must clear the CCD before the next 
exposure. Since the output register cannot be read during a row transfer (IAG 
clocking), using all the pixels, the ideal minimum read-out time is 3.685 ms. 
Placing the image in the lower left portion of the array can increase the frame rate 
since only a portion of the array needs to be read-out. After shifting the rows 
with information to the output register, the IAG clocking can increase to its 
maximum rate to clear the CCD and to commence another integration. The six 


leading dummy pixels will always add some time to each serial read-out time, but 
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the twelve trailing dark pixels could be chopped off. As an example, a 60 X 60 
image would take 2 X 60 rows X 0.667 us + 66 pixels X .1 ps X 60 rows for the 
serial read-out for a total of 476 ps. There are 66 pixels per row because of the 
six dummy pixels. 

As mentioned earlier, the integration or exposure time adds to the read- 
out time to produce the total measurement time. The total time could be reduced 
to only the integration time if the CCD were read-out while integrating the next 
set of data. This reduces the minimum total time to the longer of these two 
processes. Frame transfer devices as shown in Figure 6 have a complete extra 
set of pixels to hold the data for read-out. The interline method is not practical 
for a star image because the gaps at the transfer pixel rows reduce the exposed 
area by a factor of two. Each pixel would lose half of its photons. This is 
unacceptable for this application. Masking off half of the array is a possible 
solution. It would not work with the TC211 chip. Clocking the IAG parallel row 
shift register would produce streaking in the upper half exposed pixels while 


reading the masked pixels. 


2. Interface Card 
The Spectra Source computer interface card contains the circuitry 
needed to transfer information from the CCD array to the computer bus. It does 
an eight bit parallel transfer to an IBM compatible computer. Figure 10 is a block 
diagram of the components of this card. The CCD sensor was discussed in 


section 1. The sample-and-hold (S&H) amplifier chip converts video voltage from 
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Figure 10. Spectra Sources Lynxx PC+ interface card block diagram. 


section 1. The sample-and-hold (S&H) amplifier chip converts video voltage 
from one pixel into a voltage pulse. It is a NEC SE/NE5537 [Ref. 20]. The 
minimum acquisition time is less than 4 ps. It is controlled by pulses from the 
chip below it, a 74LS123P [Ref. 21], which is a dual retriggerable one-shot. Its 
pulse switching characteristics are in the nanosecond (ns) range. A 1200 pF 
capacitor and 12 k resistor determine the present pulse length of 5180 ns. The 
control chip AviGalaisios did not limit the speed of the S&H chip. 

The analog digital converter (ADC) chip converts the voltage to a 
positive 12 bit integer, which has a maximum value of 4095. This system has a 


Maxim MX7572 ADC chip [Ref. 22]. Its minimum conversion time is 5 ps. It, as 
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and-hold and digital conversion time is 9 ps for each pixel. For a 60 X 60 pixel 


image, the total S&H and ADC time is 32.4 ms. 


3. Computers 
Two different computers were used to investigate the dependence of 
the frame rate on computer characteristics. The first was a Compaq 386, 20 
megahertz machine. It had a 16 megabyte RAM and a 60 megabyte hard disk. 
The second was a Dell 486, 83 megahertz machine. It had a 16 megabyte RAM 


with a 350 megabyte hard disk. 


B. SOFTWARE 

The Lynxx system came with software programs designed for amateur 
astronomical observations [Ref. 17]. This thesis used a separate set of data 
acquisition program modules written for NPS by Spectra Sources [Ref. 18]. It 
was necessary to modify these modules during this thesis research. The 
software was written in the C language. Microsoft Quick C 2.5, Borland Turbo 
C++ 1.0 and Borland C++ 2.0 compilers produced executable code that were 


tested to compare their relative speeds while collecting and processing data. 
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IV RESULTS 


The first task involved writing the modulation transfer function, centroid jitter 
and display programs. After it become clear that the Lynxx PC + CCD system 
was too slow for effective centroid measurements, the second task developed 
programs to analyze the Lynxx image acquisitions times. The data acquisition 
functions were timed using different computers and code from different 
compilers. The manufacturer's maximum speed specifications for both the CCD 
camera and the interface card were compared to the actual measured times. 


After analyzing these results, replacements were recommended. 


A. PROGRAMS 

A modulation transfer function program and a jitter program were written to 
compute the atmospheric coherence length, r,. An artificial star program was 
developed to test the MTF program. The following sections describe these 


programs. 


1. Modulation Transfer Function 
This program calculated the coherence length, r,, from a stellar image 
intensity distribution. It subtracted the background from the star image, 
calculated the line spread function, found the centroid and calculated the fast 


Fourier transform. It then calculated r, using equation (9). The calculations were 
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done in the x and y dimensions separately, giving a quality check, since actual 
image turbulence is nearly always isotropic. With the image centroid data the 
telescope could operate in an auto track mode, a feature that comes with the 
Lynxx software package. The program plotted the modulation transfer function 
as a visual check. The program separates into modules for ease of future 


revisions. This program is in Appendix B. 


2. Artificial Star Program 
This program created a perfect exponential star image for the MTF 
program. This provided a convenient test of the MTF program, since the FFT of 


an exponential is also an exponential. This program is in Appendix C. 


3. Jitter 
This program calculated r, from the centroid jitter. It first removed the 
background by subtraction. After this, it computed the image centroid jitter using 
the standard deviation of the image, according to equation (11). This program 
also calculated the line spread functions and plotted them on the screen so that 


proper telescope tracking was verified. Appendix D contains this jitter program. 


B. TIME MEASUREMENTS 

Programs were written or modified to compare frame rates attainable for 
different subframe sizes. The computer clock measured the time differences. 
Since the MS-DOS system clock had a 55 ms resolution, many of the times were 


done for 100 cycles of the function of interest. All the times used were at least 
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one second long, resulting in an accuracy within 1%. An array of 5 or more trials 
ensured consistency. Appendix E is an example of the code used for measuring 


the times. 


1. Computer 

We found that the type and speed of the computer did not contribute 
significantly to the camera frame rate. Measurement time included the clearing 
of the CCD, the integration, and digitization times. With an image size of 60 X 
60, there was less than % of a frame per second difference between 386 and 
486-based computers. Figure 11 shows a comparison of frame rates verses 
subframe size for the two machines. Including calculations for the line spread 
functions and backgrounds, the frame rate improvement was about 31% frames 
per second for a 60 X 60 image with a 486 computer over that achievable with 
a 386 machine. This is shown in Figure 12. The mathematical calculations were 
timed to determine their contributions to the total frame rate. These calculations 
could have been done separately if they had slowed the frame rate, but they 


were actually performed in real time. 


2. Compiler 
It was found that the choice of compiler did not significantly influence 
the speed of the measurements. Figure 13 shows there was no difference in the 
frame rate for the code produced by three different compilers to acquire data. 


Figure 14 shows that after adding the mathematical calculations to the acquisition 
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Figure 11. Spectra Sources Lynxx PC+ frame rates using IBM PC 
386 and 486 computers. Times include clear, integration and 
digitization, using code compiled under Turbo C++. 
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Figure 12. Spectra Sources Lynxx PC+ frame rates using IBM PC 
386 and 486 computers. Times include clear, integration, 
digitization and mathematical calculations, using ccde 
compiled under Quick C. 
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COMPARISON OF COMPILERS WITHOUT CALCULATIONS 
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COMPARISON OF COMPILERS WITH CALCULATIONS 
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Figure 14. Spectra Sources Lynxx PC+ frame rates using Turbo 
C++, Borland C++ and Quick C compilers. Includes clear, 
integration, digitization and mathematical calculation times, 

using a Dell 486, 33 MHz PC. 
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time, the code produced by Quick C was about % of a frame/second slower than 


the code from Borland C++ or Turbo C++ for a 60 X 60 image size. 


3. Functions 
The times for each measurement were divided into four steps. Table 
1 shows these steps. 


TABLE 1 


The digitization and calculation times depended on the size of the image, while 












the clear and integration times were constant for each measurement. This is 
shown in Figure 15. Clear and digitization times depended on the image sensor 


and interface card. 


C. HARDWARE RATE COMPARISONS 

This section compared the maximum design speeds of the 
chips in both the camera head and interface card to those actually measured. 
A general time reference was 200 frames/second which corresponded to 5 ms 


per frame. 
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TIME VS FRAME SIZE 
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Figure 15. Spectra Sources Lynxx PC+ function times for a 
Dell 486, 33 MHz PC. Clear and integration times were 


constant at 6.45 and 5.0 ms respectively. Digitization and 
calculation times depended on the subframe size. 
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1. CCD Camera 

The measured time for the clear function was 6.45 ms. The present 
sensor, [C211, is capable of clearing the entire array in 110 ps. All rows need 
to be cleared even when reading a sub array such as a 60 X 60 image. This 
function was much slower than the sensor capability. It is not Known why the 
clear times are so slow. 

The present sensor could read a 60 X 60 image in 486 ps at maximum 
design speeds. The measured digitization time, which included the read-out 
time, sample-and-hold (S&H) and the digital conversion (ADC) times was 86.1 
ms. Since this was much larger than 436 ps, the acquisition time depended 


primarily on the S&H and ADC processes. 


2. Interface Card 
The interface card controls the S&H and ADC times. Their maximum 
speed specifications would result in processing a 60 X 60 image in 32.4 ms. This 
is a significant portion of the observed 86.1 ms digitization time. We therefore 
concluded that these two chips were the primary contributors to the frame rate 
and were not operating at their maximum speed specifications. The sensor itself 


can hypothetically operate at significantly higher speeds. 
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D. RECOMMENDATIONS 
1. CCD Camera 


To increase the frame rate requires an image sensor that uses a frame 
transfer method with separate image and storage arrays. The Texas Instruments 
TC277, a 735 X 580 pixel CCD image sensor is a proposed replacement sensor 


[Ref.19]. Figure 16 shows its set-up. The image area is 699 X 288 which is more 
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Figure 16. Recommended image sensor CCD, Texas Instruments TC277 block 
diagram. 


than big enough for our applications. The image sensing area and image 
Storage area have different gates. The image area gate (IAG) shifts the rows from 


the sensing area to the storage area one row per clock pulse. The storage area 
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gate (SAG) shifts the rows down to the serial register at one row per clock pulse. 
The three serial register gates (SRG) shift out the pixels in each row. 

The storage area gate (SAG) maximum clock rate is 3.84 MHz, which 
corresponds to 0.299 ps per row. The serial registers gates (SRG) maximum 
clock rates are 4.46 MHz, which corresponds to 0.224 us per pixel. For a 60 X 
60 image, noting that there are eleven dummy pixels per row, the read-out time 
is 60 rows X 0.299 ps + 71 pixels X 0.224 ps X 60 rows for a total of 972 us. 
This is done in parallel with the integration. Since this is much less than any 
expected integration time, the integration time and transfer time of data from the 
sensing area to the storage area become the limiting factors for frame rates. 

The image area gate (IAG) maximum clock rate is 3.34 MHz, which 
corresponds to 0.299 ps per row. The clear time would then be 288 rows X 
0.299 us or 86 us. This clear time plus the number of rows X 0.299 us plus the 
integration time will be the measurement time for each set of data. For an image 
with 60 rows, the measurement time would be only 86 + 17.9 = 104 us, plus the 
integration time. For the TC211, the minimum clear and 60 X 60 image read-out 
time is 110 + 4386 = 546 ps. The TC277 is over 5 times faster. 

It should also be noted that the number of pixels in each row of the 
TC277 can be the maximum available (735) without affecting the speed, due to 
the parallel read-out and integration times. The noise equivalent signal is only 
25 electrons, a 6-fold improvement from the TC211’s 150 electrons. Replacing 


the image sensor would require modifications to the interface board. 
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Using such an improved sensor, the bulk of the data could be 
transferred while collecting the next data set. What remains is to get the 
sampling and digital conversion times fast enough to be completed during the 


integration time. 


2. Interface Card 

The first step in the digital conversion is the sample-and-hold (S&H). 
A proposed replacement chip for this is Analog Devices AD781 [Ref. 23]. It 
provides a 700 ns conversion time. This is over 5 times faster than the present 
SE/NE5537 S&H chip's 4 ps. For a60 X 60 image, the AD781 S&H time reduces 
to 2.52 ms from 14.4 ms for the present chip. The present control chip, 
OM74LS123P, for the S&H chip does not need replacement. Changing its 
capacitor to 55 pF and resistor to 25 k would adjust its pulse width to 700 ns 
(Ref. 21]. This is the only revision necessary to replace the S&H chip. 

A proposed replacement chip for the ADC is a Maxim Max162 [Ref. 22]. 
This chip is almost identical to the present one, except that its conversion time 
is 3 ps instead of 5 ps. For a 60 X 60 image, the improved ADC time would be 
10.8 ms. The Max162 could replace the Mx7552 with little revision to the rest of 
the circuit. 

There are some faster S&H and ADC chips available if the board were 
rebuilt. Datel has a S&H chip with a time of only 25 ns and a ADC chip with a 
conversion time of 350 ns [Ref. 24]. A better option is to have the S&H and ADC 


processes combined in one chip. Of those available, the Datel ADS118 [Ref. 24] 
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has a throughput rate of 5.0 Mhz, which corresponds to only 200 ns per pixel. 
The DM74L$1283P control chip will still work by using a 60 pF capacitor with a 5 
k resistor [Ref. 21]. Besides providing fast timing, combining these two chips 
would also simplify the surrounding circuit. For a 60 X 60 image, the S&H and 


ADC time would be 720 us for the ADS118. 


E. SUMMARY 

Table 2 summarizes all of the relevant times. The components are replaced 
going from left to right across the table. They are in order of the complexity of 
replacement. Replacing the CCD sensor would be the significant step requiring 
major modifications to the interface board. There are also parallel times after the 
CCD sensor is replaced. The processes of clearing, integration and image area 
gate (IAG) transfer occur parallel to the read-out and digital conversion. The 
slowest of these times will limit the frame rate. The S&H and ADC times are the 
limiting factors. For this reason and the necessity to modify the board after 


replacing the CCD sensor, it is best to build a new system. 
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V CONCLUSIONS AND RECOMMENDATIONS 


A CCD camera can measure atmospheric turbulence by measuring the 
image spread of a point source and by measuring the centroid jitter induced by 
turbulence. The Spectra Sources Lynxx PC+ system evaluated in this thesis can 
measure the coherence length using the point source image spread technique, 
but it is much too slow for jitter measurements. Centroid motion and Greenwood 
frequency measurements need a sample rate of 200 Hz or more to avoid 
undersampling the atmospheric dynamics. Achievement of these rates requires 
a frame transfer CCD, so that parallel image exposure and frame read-out are 
possible. Data processing can be performed during an exposure if the time 
needed is sufficiently short, or after digitizing a series of exposures. The Lynxx 
PC+ system, which provides 12 frames/second, is much too slow for this. Since 
both the CCD sensor and digital conversion components need replacement, it 
is best to design and build anew system. The recommended CCD image sensor 
is a Texas Instruments TC277 [Ref. 19]. The recommended replacement for the 
sample-and-hold, and analog digital converter chips is a combined sampling 
analog-to-digital converter chip, Datel ADS118 [Ref. 24]. The atmospheric 
coherence length programs developed in this thesis should work with a new 


system. 
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The Image sensing area consists ot (65 horlzontat image lines each containing 192 ptotosensitive elements 
(pixets) Each pixetts (3.75pm (horlzontal) by 16 OOj1m (vertical) As tight enters the siticon in tlre Image sensing 
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the serial register gate can be clocked imitit all of the charge packets are moved out of the serlal register to the 
charge detection node at the amplifier Input. 


There are t2 dark pixels to the right ofthe 192 tnage ptxets on each tnage line These dark plxets are shlekded 
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Figure 1. Charge Accumutatton and Transfer Process 
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Figure 2. Timing Diagram, Noninterlace Mode 


absolute maximum ratings Over operating free-alr ternperature range (unless otherwise noted) 


Supply voltage range (ADB) (sea Noto t) .................... a Sone pas OVto t5V 
Clock voltage range for [AG, SRG, ABG terminats Soak do ee ~15Vto5V 
Operating hee-air temperature range 2.0... ccc cece eee eee e eee eteeeeees -30°C to 85°C 
SUCUGTOE ISU UIE CU RUG 1] (5s aaa 5. 5 a ~30°C to 85°C 
Lead temperature 1,6 nin (1/16 inch) from case for 10 Beco se Cee ee oe 2 260°C 


NOTE $1 All voltage vahies are with respact lo the aubatinta lorninad 
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recoinmended operating conditions 


Sipply voltage, ADB 


ame ee ee 


Sibstrale bias voltage 




















thigh bovel 
Image area gulo latormudtiato levol 
t ow tevel 
Clock vottage! Suttal rugistur gate loth 
low lovet 
ee -* High levol 
Anldstouning gale tifurmucdiute level 
1 ow luvol 
ma 1AG pe Se ee ia — 
Clock talus SRG = 7 os, Sade Teas a, i 
ae 2. ae -. 
1; Time intorval, SRG 1 to taGT 
17 Tine inforval, IAGT to toSRGtanslerputse Fo i 
13 Putso chration, [AG high 
14 Pidse dination, SAG transtor pulse igh a _ 
icc. a ee Tene wilorval, (AG! to SAG translor puso | 
i eee Tune itorval, SRG SRG transtur putse J to SIG clock piso T 
Capacitive load OUT pin Par a ae 











Opeoratiy hou ew wit lonyperalnre, Ta 





‘ieee The Bai Pees convorntion, in which the toast posthve (most noyalvo) vito by dovignalod ndndawin, s teod in fie dala shoul for clock volaye 


buvuls 
y Achustiignal ts soquinad fot oplinal petlomancy 





TEXAS wy 
INSTRUMENTS 


POST OF F1CE BUX 65530) © UAILAS TEXAS 75205 


46 


217 


TC2ti 
192 x 165-PIXEL CCD IMAGE SENSOR 


03271, JANUARY 1990 ae 


electrical characteristics over recommended operating range of supply voltage, Ta = — 10°C to 40°C 


Dynamic range (see Nole 2) 


Cherge conversion factor 

Cherge transite: efficiency (see Nolte 4) 

Signal response dolay time, 7 (826 Note § end Figure 5) 
Gamme(seeNole8) 

Output resislance 


: 14 noise (5 kHz) VTE 
Noise voltage Riendom nolse, | = 100 kHz WAAL 
Noise equivalent signal electrons 


Rejection ratio al 7 16 Miz 





From SAG to output (908 Nole 8) 
Supply current, loo 

Image aree gale 
Capacitance Sorinl ragisier gate 





PARAMETER 
Anliblooming disabled (see Nofe 3) 
Anlibloorning enabled 












From ADB lo outprd (soe Note 7) 








Antiblooming gate 


Tan lypical values ere at Ta = 25°C 








NOTES: 2) Dynamic tango is — 20 tinos tha logarittiin of tho monn noise signal divided by the salutation output signal 
3 For this test, the entdldooming gale mirst bo hinend al tho lnlonnedialn boval 
4° Charge transler efficiancy bs one minis the charge toss per transtear in tho ontpul register Tha lost le porformed hi the dark uelng an 
eloctrical Inpirt signal 
§ Signal msporsa doloy Inne le the tine hotwoon Ihe falling odgn of ha SAG clock pulse arul the outpd signal valid slato 
6 Garona {y} Is fre vate of the exponont in Wha equation below fortwo points on fre Buest portion of Ile lranslor hunction curve {this 
vahis represents points naer salrration) 
Expostme (2) 7 | Ortterdt signal (2) 
Exposine (1) Oudpet signal (1) 
7 ADB rejection ralo fs - 20 Umes Ihe logaidtien of tho ac anphlades al Ihe output divided by fro ac ainphhide al ADB 
B SAG tajoction ratio le - 20 timros Ihe logarithe of the ac angilucle al Iho oitpil drvided by tho ac anphlide at SAG 
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2-16 FOS! OF FICE BOX 855903 © DAILAS, TEXAS 7526S 


47 


TC2t1 
192 x 165-PIXEL CCD IMAGE SENSOR 


03271, JANUIARY 1990 





optical characteristics, Ta = 25°C (untess otherwise noted) 


PARAMETER 
No IA fitter 
With IR ther 


Seusiwity (see Note 9) Muasinud al Viyy (suo Noto 10) 


Antiblooming disablod 
Antiblooming onablud 
Strabe 

Shuttorod kot 


Saturation signal (see Note 11) 


Blooning overload rato (soe Note ¢t 2) 


Output signal nonundormity (1/2 seturation) (see Note 13) 
Image atoa welt capacity 


Derk cntrent Ta=20°C ; 


Dark signal (see Note 14) 

Dark signal nonuniformity tor ente hold (see Nolte 15) 

sagt eS Horizontal 
Vortic al 


NOTES: 9 Sonsfivity le tneasured et an inlegt ation tino of 16 667 nts end © soinco tumporatino of 2856 K ACM 500 titor is used 
10 Vi, is tho outrad vottuge thal ropresonts thy Mtnostokl of oporation of ontiddooming Vey ~ t/2 salaration signel 
11 Saturetion Is the condition In which furthot bicreuso In exposito dues not fuad to futhor incruaso in ordynd signal 
12 Blooinintg ovurload tubo is the ruho of tlooming expositro to eaturalion oxposnte 
13 Gutpat signal nonnmiorniity ts the ratio of tho maxtnuin ys: ol fo pizul chorunce hi outpet signal to tho moan ouad signal for exposriie 
8 pistud tu give 1/2 thy satinubon outpul signal 
14 Dark signal fovel rs mnoasired fom the dummy piruls 
15 Dark stgnol nonmilormiy ts tho maximum pizul to pleul dituronce ite dak comblion 


{Aodulution transfer hincton 





PARAMETER MEASUREMENT INFORMATION 


Vinicry min 





Figure 3. Typical Clock Waveform for [AG and ABG 


SLEW RATE BETWEEN 10% and 00% = 70 to 120 Vis. t, = 150 ne, t) - 90 ne 
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PARAMETER MEASUREMENT INFORMATION — 


VHIGH min 





SLEW RATE BETWEEN t0% and 00% = 300 Vis, t, = tge f5 ns 


Figure 4. Typical Clock Waveform for SRG 
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Figure 5. SRG and CCD Output Waveforms 
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MTF — Modiutation Traenster Function 


TYPICAL CHARACTERISTICS 


HORIZONTAL MODULATION 
TRANSFER FUNCTION 
(BARS PARALLEL TO SERIAL REGISTER) 


1 | ee a ——_ 
—— 


086 





02 h = 400 to 700 nm Monuchromatic tighe 


MTF —- Mocuiston Transter Function 
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Figure 6 


NOISE SPECTRUM OF OUTPUT AMPLIFIER 
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VERTICAL. MODULATION 
TRANSFER FUNCTION 
(BARS PERPENDICULAR TO SERIAL REGISTER) 
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Figure 7 


CCD SPECTRAL RESPONSIVITY 
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TYPICAL APPLICATION DATA 
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Figure 10. Typical Application Circult Dlagram 
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MECHANICAL DATA 


e } ; 
Thte package for the TC21t consists of a ceramic base, glass window, and a 6 lead frame. The glass window 


is Sealed to the package by an epoxy adhesive The package leads are configiuredin a dual in (ine organization 
and {il into mounting holes with 2,54 min (0 1 inch) center-lo cenler spacings 


_..p] _ 7,64 (0 297) 
7,14 (0 261) 








(0 012 A 7,62 {0 300] --- 
fo03) 


whe Saeore 
2,54 {0 100} 0,36 (0 Ot5) 
30 


NOTES A Dunenstons ore ty mifinetors and parontiiete ally mi hus Singh dunerstons ure nonin 
B Theconer of the package und the contur of Ku nage ios Bre nob conn thor 


C Thodistenco trom the top of the ghtss ly the imaye sensor sttface ts lypr tty Core (0 040 tech) Thy ghise is lypicully 0 020 iict thick 
andhas an indor of teftaction of | 52 
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APPENDIX B MTF PROGRAM 


eMF.C find ro and f from the Isf using FFT to the MTF*******/ 
#define TITLE "Measure ro, plot the MTF" 
#define AUTHOR "By W. J. Rall" 


ee Fee REAL DOK Cee an eee SR ae See ee ARR 


/*#include <alloc.h> */ 
#include <lynxx.hn> 


#define True 1 
#define False O 
#define false O 


#include <stdio.h> 
#include <stdlib.h> 
#include <graph.h> 
#include <math.h> 


#include <sys\types.h> 


/*int frame_ptr[FRAME_BYTES/2]:*/ 
int size=60, xstart=0O, ystart=O; 
float Isfx[165],lsfmax=0.0,Isfy[192]; 


int sign=1, m=8, mm=2, mm2:/*mm is # of points in FFT, mm2=mm/2***/ 
float mtfx[1024], mtfmax=0.0,mtfy[1024], re[1024],im[1024],px[512], py[512]; 
float rocal = 1.0,rox,roy; 


void init(): 
void subexpose(long tm,int shutter); 
void expose(long tm); 
void long_delay(long tm) ; 
void calc Isf(unsigned int far *sub_ frame); 
void background(); 
void fft(int m, int sign, float re[], float im[]}); 
void calc_ro(); 
void init_mtf(); 
void plot_mttf(): 
/*to see Isf's, use these two functions instead of mtf ones*/ 


53 


void init_graph(); 
void plot_graph(); 


[B------------------------ main routine----------------------- 7 
void main() 
{ . . 
Int I: 
for (i=0; i<m-1; 1++) mm = 2*%mm; 
init_lynxx(); 
init(); 


cooler_off(); 
_setvideomode( _TEXTC80); 


} 
[*-------------------- Initialize input info------------------- =) 
void init() 
{ 
int cool = O,shutter=O0; 
long tm; 


printf("Expose time in ms.:"); 
scanf("%ld", &tm); 
printf("Cooler O=OFF 1=ON:"); 
scanf("%d",&cool); 
if (cool) 
cooler on(); 
printf("Use shutter ? (1=Yes/O=No)'); 
scanf("%d", &shutter); 
printf("Subframe size? "); 
scanf("%d", &size); 


/* printf(\nHit any key to quit:");*/ 


subexpose(tm,shutter) ; 


} 

[* ------------------------ subExpose CCD------------------------ 
void subexpose(long tm,int shutter) 

{ 


unsigned int *sub_ frame; 
int sub_frame_bytes; 


float px0,py0; 
Int 1,},X,Y; 
for (i=0; i<m-1; i++) mm = 2*%mm; 


mm2= (int) ((float) (mm)/2.0): 


sub frame_bytes = size*size * sizeof(unsignedé int); 
sub frame = (unsigned int*)malloc(sub_frame_bytes); 


if (!shutter) open_shutter(); 


while ( !kbhit() ) 


{ 
if (shutter) 


expose(tm); 


} 


else 


Seca,  umnor USING SHULer © Fhe *** tse enn ann sx | 
long delay(tm); 

} 

xX param = (unsigned char) xstart; 

y_param = (unsigned char) ystart; 

size param = (unsigned char) size: 

ptr param = sub frame; 

digitize sub frame(); 

calc_Isf(sub_ frame); 

background(); 


for (x=xstart+size-1; x>=xstart; x--) re[x] = Isfx[x]; 
fft(m, sign, re, im); 
for (x=0; x<mm; x++) mtfx[x] = re[x]; 


aAManne xo OOWer SPECIFUITN® "== - 5" =" "** **/ 
px0O = re[O]*re{O] + im{[O]*im[O}; 
for (j=0; j<mma2; j++) 
{ 
if(mtfxfj] != 0.0 ) if(px0 != 0.0) 
px{j]=sqrt((re[j]*reY] +im[j] *im[j])/px0)/mttx{j]; 


} 
/*reinitialize the FFT arguements before finding the Y FFI*/ 
for (i=O; 1< 1024; i++) 


{ 
re[i] = 0.0; 
imf[i] = 0.0; 
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for (y=ystart+ size-1; y>=ystart; y--) re[y] = Isfy[y]; 


fft(m, sign, re, im); 
for (y=0; y<mm; y++) mtfy[y] = rely]; 


/*find the y power spectnimit** 273242 eas, 
pyO = re(0]*re[O] + im[O]*im[0]; 
for (j=0; j<mm2; j++) 
{ 
if(mtfy[j] != 0.0) if(pyO != 0.0) 


pyU]=sart((refj]*re[j] +im[j]*im[j])/py)/mtty [j); 


calc_ro(); 

init_mtf(); 

plot_mitf(); 

/*init_graph(); 

plot graph(); ********* use these to see Isf’s 


free(sub_ frame); 
close_shutter(); 


} 
| * ------------------------ Expose CCD------------------------ =) 
void expose(long tm) 
{ 
CIKECEe 


open_shutter(); 
long _delay(tm): 
close_shutter(); 


} 

[*------------------------ long delay------------------------ | 
void long delay(long tm) 

{ 


delay param = 1000: 
while (tm > 1000) 


{ 


delay_ms(); 
tm = tm -1000; 
} 


delay param = (unsigned short)tm; 
delay _ms(); 


void calc_Isf(unsigned int far *sub frame) 


{ 
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unsigned int pixel; 

Int x,y; 

float Isfkmax=0.0,lsfymax= 0.0; 

for (x=0; x<165; x++) Isfx[x]=0.0; 
for (y=0; y<165; y+ +) Isfy[y]=0.0; 
for (x=size-1+xstart; x>=xstart; x--) 


{ 

for (y=size-1+ystart; y>=ystart; y--) 

{ 
pixel=*(sub_frame+ +); 
Isfx[x] = Isfx[x] + (float)pixel; 
Isfy[y] = Isfy[y] + (float)pixel; 
if(Isfx[x] > Isfxmax) Isfxmax=Isfx[x]; 
if(Isfy[y] > Isfymax) Isfymax=Isfy[y]; 

} 


Isfmax=Isfxmax; 
if(isfxmax < Isfymax) Isfmax = Isfymax; 


void background() 
{ 
int x,y,edge; 
float sumx=0.0,sumy=0.0,ave x1,ave_x2,ave yl,ave y2,aveback; 
edge=(int)(size/10.0); /*this uses 20% to find background*/ 
/*-----Find the average background--------------- ai | 
for (x=xstart; x<xstart+edge; x++) sumx=sumx-+lsfx[x]; 
ave _x1=sumx/edge; sumx=0.0; 
for (x=xstart+ size-edge-1; x<xstart+size; x++) sumx=sumx+lIsfx[x]; 
ave _x2=sumx/edge; 
for (y=ystart; y<ystart+ edge; y++) sumy=sumy +lsfy[y]; 
ave y1=sumy/edge; sumy=0.0; 
for (y=ystart+size-edge-1; y<ystart+size; y+ +) sumy=sumy +lsfy[y]; 
ave _y2=sumy/edge; 
aveback =((ave x1+ave_ x2)+(ave y1+ave_y2))/4.0; 
/*----subtract off the background------- ad 
Isfmax = Isfmax - aveback: ; 
for (x=xstart+size-1; x>=xstart; x--) 
Isfx[x] =Isfx[x]-((x-xstart) *(ave_x2-ave_x1)/(size-edge)+ave_x1); 
for (y=ystart+ size-1; y>=ystart; y--) 
Isfy[y]=Isfy[y]-((y-ystart)*(ave_y2-ave_ y1)/(size-edge)+ave_y1); 
/* printf(’\nThe average background is %f\n",aveback);*/ 


of 


} 


[RR RRR AAANAAD Dee oceans the FP] MARRS SHH RS AR ee 
void fft(int m, int sign, float re[], float im{]) 
{ 

int n, j, Jj1, ndiv2, n1, n2, k, i, ip, npts; 

int le, leO, le1, I; 

float pts, t, ure, uim, wre, wim, tre, tim; 

double ang, pi; 


pi = 4.0*atan(1.0); 
is Aik ills td (double)m) + 1.0e-10): 


n2=n- 2 
ndiv2 = n/(int)2; 


for(i=O; i< =n2; i+ +) 


{ 
if(i < j1) 
{ 


t = re[j1]; /* you might do better here with pointer operations 
re(j1] = re[i]; 
re[i] = t; 
t = im{j1]; 
Iimfj1] = imf{i]; 
Im{i] = t; 
} 
k = ndive: 
while(k < j) 
{ 


] -= k; 
k /= 2: 


} 


le = 1; 
for(l=1; I<=m; I++) 


{ 
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ure = (float)1; 

uim = (float)O; 

ang = pi/(double)le; 
wre = (float)cos(ang); 
wim = (float)sin(ang); 


leO=le; 
le1=le-1: 
le +=le; 


for(j=O; j<=le1; j++) 


{ 
for(i=j; i<=n1; 1+ =Ile) 
{ 
ip = | + led; 
tre = re[ip]*ure - im[ip]*uim; 
tim = re[ip]*uim + im[ip]*ure; 
re[ip] = re[i] - tre: 
im[ip] = im{[i] - tim; 
re(i] += tre; 
im[i] += tim; 
} 
t = ure*wre - uim*wim: 
uim = ure*wim + ulm*wre; 
ure = ft; 
} 
F 
if(sign<0O) 
pts = 1/(float)n; 
for(i=O; i<n; i++) 
if 
re(i] *= pts; 
im{i] *= pts; 
} 


ee AR KAR AA tO Olin aie ne OR Re AAA K KK / 


void calc_ro() 


int j; 
float z,z0ld,pj,slope; 


if (px[1] > 1.0) if (px[1] <= 0.0) printf("\nX DATA POINT BAD"): 


3S 


else 
{ 
2. = log (-log(px[1])); 


{ 
Pp} = pxij]; 
zold = Z; 
if (pj<1.0) 
z = log(-log(pj)); 
slope = (z-zold)/log( (float) (j)/(float) (j-1)); 
if (2 > 0.0) goto ro; 
} 


} 


ro: rox = rocal*exp(-zold/slope + log(j-1)); 


if (py[1] > 1.0) if (py[1] <= 0.0) printf("\nY DATA POINT BAD"): 
else 


{ 
Z = log (-log(py[1])); 
for j=2; j<mm2; j++) 


{ e e 
Pp) = pylj]: 
zold = 2; 
if (pj <1.0) 
z = log(-log(p))); 
slope = (z-zold)/log( (float) (j)/(float) (j-1)); 
if (z > 0.0) goto ro2:; 
} 
} 
} 
ro2: roy = rocal*exp(-zold/slope + log(j-1)); 
} 
[* ------------------ initialize for plotting--MTF’s------------------- */ 
void init_mtf() 
int TRUE=1; 
int x,y; 


for (x=0; x<mm; x++) 
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if(mtfx[x] > mtfmax) mtfmax=mtfx[x]; 
for (y=0; y<mm; y+ +) 
if(mtfy[y] > mtfmax) mtfmax=mitfy[y]; 
_setvideomode( VRES16COLOR); 
_setviewport(0,39,639,439); 
_setwindow(TRUE,0,0,mm+mm,mtfmax); 


_setcolor(15); 
_moveto_w<(0.,0.); 


void plot_mtf() 
{ | 
int x,y; 
for (x=mm2; x<mm; x++) 
_lineto_w((double) (x-mm2), (double) mtfx[x]); 


for (x=0; x<mm2; x++) 


_lineto_w((double) (x+mm_2), (double) mtfx[x]); 
} 


for (y=mm2; y<mm; y+ +) 
_lineto_w((double)(mm2+ y),(double)mtfy[y]); 


for (y=0; y<mm2; y++) 


{ 
} 


_settextposition(1,1); 
printf(’%s\n %s", TITLE, DATE_); 


_lineto_w((double)(mm+ mm2 + y),(double)mtfy[y]); 


_settextposition (4,1); 

printf("rox = %f",rox); 
_settextposition (4,38); 
printf("roy = %f" roy); 


_settextposition (40, 30); 
_outtext("Hit any key to exit:"); 
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/* getch(); there is also a getch in subexpose, only use one*/ 


void init_graph() 


int (RUE=f 
_setvideomode( VRES16COLOR): 
_setviewport(0,39,639,439); 
_setwindow(TRUE,0,0,size+ size, lsfmax): 
_setcolor(15); 
_moveto_w(0.,0.); 


void plot_graph() 
{ 
Int X,Y; 


for (x=xstart; x<xstart+size; x++) 


{ 
_lineto_w((double) (x-xstart), (double)|Isfx[x]); 


} 


for (y=ystart; y<ystart+size; y++) 
{ 
_lineto_w((double)(size+ y-ystart), (double)|sfy [y]); 
} 


_settextposition(1,1); 
printf("%s\n %s\n\n", TITLE, DATE_): 
_settextposition (40, 30); 
_outtext("Hit any key to exit:"); 
/* getch(); — there is also a getch In subexpose, only use one*/ 
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APPENDIX C ARTIFICIAL STAR PROGRAM 


/*Uses an artificial star to test the MTF.c program*has Isf plot also**/ 
#define TITLE "Plot an artificial star MTF" 


#define AUTHOR "By W. J. Rall" 


eer = CUNCK Ce ea ee ee RN KK KKK | 


#include <stdio.h> 
#include <stdlib.h> 
#include <graph.h> 
#include <math.h> 


#include <sys\types.h> 


int xsize=165, ysize=192, xstart=0, ystart=0; 
float Isfx[165],lsfmax=0.0,Isfy[192]; 


int sign=1, mM=8, mm=2, mm2; 
float mtfx[1024],mtfmax=0.0,mtfy[1024], re[1024],im[1024],px[512],py[512]; 
float rocal = 1.0,rox,roy; 


void make _star(); 

void background(); 

void fft(int m, int sign, float re[], float im[]); 

void calc _ro(); 

void init_mtt(); 

void plot_mttf(); 

/*to see Isf's use these two functions instead of mtf ones*/ 
void init_graph(); 

void plot_graph(); 


void main() 
{ 
float px0,py0; 
Int 1,},X,Y; 
for (i=0; i<m-1; i++) mm = 2*mm; 
mm2= (int) ((float)(mm)/2.0); 
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make_star(); 
background(); 


for (x=xstart+xsize-1; x>=xstart; x--) re[x] = Isfx[x]; 
fft(m, sign, re, Im); 
for (x=0; x<mm; x++) mtfx[x] = re[x]; 


/*tind the x power spectrulmit>* "++" #45 anos) 
pxO = re[O]*re[O] + im[O]*im[O]; 
for (j=0; j<mm2; j++) 


if(mttx[j] != 0.0 ) if(pxO != 0.0) 
px[j]=sqrt((re[j]*re[j] +im[j] *im{j])/px0)/mttx{]; 


/*reinitialize the FFT arguements before finding the Y FFIT*/ 
for (i=O; i< 1024; i++) 


refi] = 0.0; 
im{i] = 0.0; 
} 


for (y=ystart+ysize-1; y>=ystart; y--) re[y] = Isfy[y]: 
fft(m, sign, re, im); 
for (y=0; y<mm; y++) mitfy[y] = rely]; 


/*find the y power spectrum*************%*%/ 
pyO = re[O]*re[O] + im[O]*im[0}; 
for (j=0; j<mm2; j++) 


if(mtty[j] != 0.0) if(pyO != 0.0) 
pyUj]=sqrt((re[j]*re[j] +im[j]*im[])/py0)/mittfy Yj]; 


calc_ro(); 

iInit_mtf(); 

plot_mtf(); 

/* — init_graph(); 
plot_graph(); a 


_setvideomode( _TEXTC80); 


void make_star() 


64 


int x,y; 
float pixel,|sfxmax=0.0,Isfymax=0.0; 
for (x=xstart+xsize-1; x>=xstart; x--) 


{ 
for (y=ystart+ysize-1; y>=ystart; y--) 


pixel =4095*exp(-0. 1 * (float) (((x-82) *(x-82) + (y-96) *(y-96)))); 
Istx(x] = Isfx[x] + pixel; 
Isfy[y] = Isfy[y] + pixel: 
if(Isfx[x] > Isfxmax) Isfxmax=Isfx[x]; 
if(Isfy[y] > Isfymax) Isfymax=Isfy[y]; 
} 


Isfmax=Isfxmax: 
if(Istxmax < Isfymax) Isfmax = Isfymax; 


} 

[*------------- oubtract off the background---------------------- | 
void background() 

{ 


int x,y,edge; 
float sumx=0.0,sumy=0.0,ave x1,ave x2,ave yl,ave_ y2,aveback; 
edge= (int) (ysize/10.0); /*this uses 20% to find background*/ 
/*-----Find the average background--------------- i 
for (x=xstart; x<xstart+edge; x++) sumx=sumx+lIsfx[x]; 
ave x1=sumx/edge; sumx=0.0; 
for (x=xstart+ xsize-edge-1; x<xstart+xsize; x+ +) sumx=sumx-+Isfx[x]; 
ave x2=sumx/edge; 
for (y=ystart; y<ystart+edge; y++) sumy=sumy +lsfy[y]; 
ave y1=sumy/edge; sumy=0.0; 
for (y=ystart+ysize-edge-1; y<ystart+ysize; yt) 
sumy=sumy +lsfy [y]; 
ave _y2=sumy/edge; 
aveback =((ave_x1+ave x2)+(ave y1+ave_y2))/4.0; 
/*----subtract off the background------- ah 
Isfmax = Isfmax - aveback: 
for (x=xstart+xsize-1; x>=xstart; x--) 
Isfx[x] =Isfx[x]-((x-xstart)*(ave_x2-ave_x1)/(xsize-edge) +ave_x1); 
for (y=ystart+ysize-1; y>=ystart; y--) 
Isfy [y]=Isfy[y]-((y-ystart) *(ave_y2-ave_y1)/(ysize-edge)+ave_y1); 
/* printf("\nThe average background is %f\n",aveback);*/ 


} 
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[REREASESS Tee ee) the FR TARE RRR RAR AR ERA AE eae 


void fft(int m, int sign, float re[], float im[]) 


{ 
Int n, j, j1, ndiv2, n1, n2, k, |, ip, npts; 
int le, leO, le1, I; 
float pts, t, ure, ulm, wre, wim, tre, tim; 
double ang, pi; 
pi = 4.0*atan(1.0); 
n = (int)(pow(2.0, (double)m)+1.0e-10); 
j= 1, 
J =o 
ni=n-1: 
n2eé=n-2: 
ndiv2 = n/(int)2: 
for(i=O; i< =n2; i++) 
{ 
ieee) 
{ 
t = re[j1]; /* you might do better here with pointer operations 
*/ 
re{j1] = refi]: 
refi] = t; 
t = imfj1]; 
Im{j1] = imfi); 
imfi] = t; 
} 
= ndiv2: 
while(k < }) 
{ 
j -= k; 
K /= 2; 
} 
J or= be 
p= Puy 
} 
le = 1; 
for(l=1; |< =m; 1++) 
{ 


ure = (float)1; 
uim = (float)O; 
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ang = pi/(double)le; 
wre = (float)cos(ang); 
wim = (float)sin(ang); 


leQ=le; 
le1=le-1-: 
le +=le: 


for(j=0; j<=le1; j++) 


{ 

for(i=j; i<=n1; i+=le) 
ip = 1 + leO; 
tre = re[ip]*ure - im[ip]*uim; 
tim = re[ip]*uim + im[ip]*ure; 
re[ip] = re[i] - tre; 
im[ip] = imf{i] - tim; 
re[i] += tre: 
im[i] += tim; 

} 

t = ure*wre - uim*wim; 

uim = ure*wim + uim*wre; 

ure = t: 

} 
} 
if(sign<Q) 

pts = 1/(float)n; 
for(i=0; i<n; i++) 
re[i] *= pts; 
im{i] *= pts; 

} 


ieee 6 Re * **X**X* COlCculate iOus aeiia ee te RARER ERK / 
void calc _ro() 


Int J; 
float z,zold,pj,slope; 


if (ox[1] > 1.0) if (px[1] <= 0.0) printf("\nX DATA POINT BAD"): 
else 


{ 
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z = log (-log(px[1])): 
for (j=2; j<mm2; j++) 


Pp) = Pxil; 
Zo\en— Zs 
if (pj< 1.0) 


Z = log(-log(pj)); 
slope = (z-zold)/log( (float) (j)/(float) (j-1)); 
if (z > 0.0) goto ro; 


} 
} 
ro: rox = rocal*exp(-zold/slope + log(j-1)): 
if (py[1] > 1.0) if (py[1] <= 0.0) printf(‘\nY DATA POINT BAD‘); 
else 
{ 
z = log (-log(py[1])); 
for (j=2; j<mme2; j++) 


{ 
P) = pyUI 
zold = z 
if (pj <1.0) 
Z = log(-log(p))); 
slope = (z-zold)/log((float) (j)/(float) j-1)); 
If (Zz > 0.0) goto ro2; 
} 
} 
} 
ro2: roy = rocal*exp(-zold/slope + log(j-1)): 
j 
[ * ------------2----- initialize for plotting--MTF’s-------------- a 
void init_mtf() 
{ 
int TRUE=1: 
int x,y; 


for (x=0; x<mm,; x++) 
if(mtfx[x] > mtfmax) mtfmax=mtfx[x]; 
for (y=0; y<mm; y+ +) 
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if(mtfy[y] > mtfmax) mtfmax=mtty[y]: 


setvideomode( VRES16COLOR); 

setviewport(0, 39,639,439); 
_setwindow(TRUE,0,0,mm+mm,mttmax); 

setcolor(15); 

moveto_w(0.,0.); 


void plot_mttf() 


{ 


int x,y; 
for (x=mm2; x<mm; x++) 


_lineto_w((double) (x-mm_2), (double) mtfx[x]); 


for (x=0; x<mm2; x++) 


{ 
} 


for (y=mm2; y<mm; y++) 


_lineto_w((double) (x+mm_2), (double)mttx[x]); 


_lineto_w((double)(mm2+ y),(double)mtfy[y]); 
for (y=0; y<mm2; y++) 


_lineto_w((double)(mm+ mm2 + y),(double)mtfy[y]); 
j 


_settextposition (1,1); 
printf("%s\n Ys", TITLE, DATE_); 


_settextposition (4, 1); 

printt("rox = %f",rox); 
_settextposition (4,38); 
printf("roy = %f" roy): 


settextposition (40,30): 


outtext("Hit any key to exit:"); 
getch(); 
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void init_graph() 


int TRUE= 
_setvideomode( VRES16COLOR); 
_setviewport(0,39,639, 439); 
_setwindow(TRUE,0,0,xsize+ysize,|sfmax): 
_setcolor(15); 
_moveto_w(0.,0.); 


void plot_graph() 
{ 
int x,y; 


for (x=xstart; x<xstart+xsize; x++) 


{ 
_lineto_w((double) (x-xstart), (double)lsfx[x]); 
} 
for (y=ystart; y<ystart+ysize; y+ +) 
{ 
_lineto_w((double)(xsize+ y-ystart), (double)lsfy [y ]):; 
I 


_settextposition(1,1); 

printf("%s\n %s\n\n", TITLE, DATE_); 
_settextposition (40,30); 
_outtext("Hit any key to exit:"); 

getch(), 
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APPENDIX D JITTER PROGRAM 


/*jitter.c**********finds the cetroid,jitter, ro and f, plots Isf******/ 
#define TITLE "Jitter calculations for ro and f, plots Isf's" 


#define AUTHOR "By W. J. Rall" 

pas eee OiCK Cras eer ae AERRAAAAA ER Ce KERR RER EAR EE / 
/*#include <alloc.n> */ 
#include <lynxx.n> 


#define True 1 
#define False O 
#define false O 


#include <stdio.h> 
#include <stdlib.h> 
#include <graph.h> 
#include <math.h> 
#include <time.h> 


#include <sys\types.h> 
/*int frame_ptr[FRAME BYTES/2];*/ 


int size=160, xstart=0, ystart=O; 
float Isfx[165],lsfmax=0.0,Isfy[192]: 


float xcsum=0.0, xc2sum=0.0, ycsum=0.0, yc2sum=0.0; /*for jitter calcs*/ 


intn=1; /*number of centroids between each jitter calc*/ 
float dia=1.0, focal=1000.0, k=12.57E+6:/*wave number, 500nm*/ 
float xjit,yjit, rocal = 1.0, rox,roy,fx,fy; 


void init():; 
void subexpose(long tm,int shutter): 
void expose(long tm); 
void long delay(long tm) ; 
void calc_Isf(unsigned int far *sub_ frame); 
void background(); 
void centroid(); 
void jit(): 
void jit_rof(); 
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void init_graph(); 
void plot_graph(); 


void main() 


{ 
init_lynxx(Q); 
init(); 
cooler_off(); 
_setvideomode( _TEXTC80); 


void init() 

{ 
int cool = O,shutter=0; 
long tm; 


printf("Expose time in ms.:"); 
scanf("%ld",&tm); 


printf("Cooler O=OFF 1=ON:"); 
scanf("%d",&cool): 
if (cool) cooler on(); 


printf('Use shutter ? (1=Yes/O=No)'); 
scanf("%d", &shutter); 


printf("Enter size of subframe in pixels:’); 
scanf("%d" &size): 


printf("\nHit any key to quit:"); 


subexpose(tm, shutter); 


void subexpose(long tm,int shutter) 


{ 


unsigned int *sub_ frame; 
int sub_frame_bytes, |; 
double deltime,time1 ,time2: 


UZ 


sub_frame_bytes = size*size * sizeof(unsigned int); 
sub_frame = (unsigned int*)malloc(sub_frame_bytes); 


if (!shutter) open_shutter(); 


while ( !kbhit() ) 


{ 


for(i=O; i<n; i++) 


{ 
if (Shutter) expose(tm); 
else 
circcd(); /*if not using shutter*/ 
long_delay(tm); 
X param = (unsigned char) xstart; 
y param = (unsigned char) ystart; 
size param = (unsigned char) size; 
ptr param = sub frame; 
digitize sub frame(); 
calc_Isf(sub_ frame); 
background(): 
centroid(); 
} 
jit(); 
jit_rof(); 
init_graph() 


plot g raph(): 


free(sub_ frame); 
close_shutter(); 


} 
| *------------------------ Expose CCD------------------------ of | 
void expose(long tm) 
{ 
cirecd(); 


open_shutter(); 
long delay (tm); 
close_shutter(); 
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void long delay(long tm) 


{ 
delay param = 1000; 


while (tm > 1000) 


{ 
delay _ms(); 
tm = tm -1000; 
} 
delay param = (unsigned short)tm; 
delay _ms(); 
I 
¥en nen wenn nn nen nnn -ne calculate the line spread functions----------- oy 
void calc _Isf(unsigned int far *sub_ frame) 
{ 
unsigned int pixel; 
int x,y; 


float Isfkmax=0.0,lsfymax=0.0; 
for (x=0; x<165; x++) Isfx[x]=0.0; 
for (y=0; y<165; y+ +) Isfy[y]=0.0; 
for (x=size-1+xstart; x>=xstart; x--) 
{ 
for (y=size-1+ystart; y>=ystart; y--) 
{ 
pixel=*(sub_frame+ +); 

Isfx[x] = Isfx[x] + (float) pixel: 
Isfy[y] = Isfy[y] + (float)pixel; 
if(Isfx[x] > Isfxmax) Isfxmax=Ilsfx[x]; 
if(Isfy[y] > Isfymax) Isfymax=lsfy[y]; 

} 


Isfmax=Isfxmax; 
if(isfkmax < Isfymax) Isfmax = Isfymax; 


void background() 


{ 
int x,y,edge; 
float sumx=0.0,sumy=0.0,ave_x1,ave_x2,ave_y1,ave_y2,aveback; 
edge=(int)(size/10.0); /*this uses 20% to find background*/ 
/*-----Find the average background--------------- “i 
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for (x=xstart; x<xstart+ edge; x+ +) sumx=sumx+lsfx[x]; 
ave x1=sumx/edge; sumx=0.0; 
for (x=xstart+size-edge-1; x<xstart+size; x++) sumx=sumx-+lsfx[x]; 
ave x2=sumyx/edge; 
for (y=ystart; y<ystart+ edge; y+ +) sumy=sumy +lsfy[y]; 
ave yi=sumy/edge; sumy=0.0; 
for (y=ystart+ size-edge-1; y<ystart+size; y+ +) sumy=sumy +lsfy[y]; 
ave _y2=sumy/edge; 
aveback =((ave x1+ave_ x2)+(ave _y1l+ave_y2))/4.0; 
/*----Subtract off the background------- a 
Isfmax = Isfmax - aveback; 
for (x=xstart+size-1; x>=xstart; x--) 
Isfx[x] =Isfx[x]-((x-xstart)*(ave_x2-ave_x1)/(size-edge)+ave x1); 
for (y=ystart+size-1; y>=ystart; y--) 
Isfy [y]=Isfy[y]-((y-ystart)*(ave_y2-ave_y1)/(size-edge)+ave_y1); 
/* printf("\nThe average background Is %f\n",aveback);*/ 


void centroid() 


int x,y; 
float xcent,ycent,sum|Ix=0.0,sumly =0.0,sumxlx=0.0,sumyly =0.0; 


for (x=xstart+ size-1; x>=xstart; x--) 
for (y=ystart+size-1; y>=ystart; y--) 


Ssumxilx = sumxlx+ Isfx[x]* (float) (x); 
Sumlx = sumIx + Isfx[x]; 

sSumyly = sumyly+ Isfy[y]*(float)(y); 
sumly = sumly + Isfy[y]; 

ii 

} 


xcent=sumxIx/sumlx; ycent=sumyly/sumly; 

/*printf("\n Centroid is xcent=%f \n ycent=%f",xcent,ycent) ;*/ 
acum centroias for jitter cales**=*************/ 

xcsum = xcsum + xcent; xc2sum = xc2sum + xcent*xcent: 
ycsum = ycsum + ycent; yc2sum = yc2sum + ycent*ycent: 


} 


(| Rl lala ee ee ee eAEX OAC late the iicuma: °° noe eee x / 


void jit() 
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xjit = sqrt((xc2sum-xcsum)/(float)(n)); 

yjit = sqrt((yc2sum-ycsum)/(float)(n)); 

xcsum=0.0; xc2sum=0.0; ycsum=0.0; yc2sum=0.0; 

printf("\n Jitter is xjitter=%f  yjitter=%f" xjit, yjit):; | 
| 


[RRRAERREARAAA AEE AT CAC ate ro and f *from the jitterA AAs Se eae , 
void jit_rof() | 





{ 
float a,b; 
: 
a = xjit * k / focal: | 
b = 1.418 * pow(dia,0.333) * pow(a,2); 
rox = rocal / pow(b,0.6); 
a = yjit * k / focal: 
b = 1.418 * pow(dia,0.333) * pow(a,2); 
roy = rocal / pow(b,0.6); 
*X = xijit; 
fy = yjit*roy; 
} 
| Soceneennneeeeeene initialize for plotting--------------------- */ 


void init_graph() 
{ 


int TRUE=1; 
_setvideomode( VRES16COLOR); 
_setviewport(0,39, 639,439); 
_setwindow(TRUE,0,0,size+ size, lsfmax); 
_setcolor(15); 
_moveto_w(0.,0.); 


void plot_graph() 
i 
int x,y; 
for (x=xstart; x<xstart+size; x++) 


{ 
_lineto_w((double) (x-xstart) , (double)|sfx[x]); 
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} 


for (y=ystart; y<ystart+size; y++) 
_lineto_w((double)(size+ y-ystart), (double)lsfy[y]); 


} 
_settextposition (1,1); 
printf("%s\n %s\n\n", TITLE, DATE_ ): 


_settextposition (4,1); 

printf("rox = %g\nfxk = %g",rox,fx); 
_settextposition (4,38); 

printf("roy = %g",roy); 
_settextposition (5,38); 

printt("ty = %g",fy); 
_settextposition(40, 30); 


_outtext ("Hit any key to exit:"); 
/* getch(); subexpose also has one, only use one*/ 


} 


ae 


APPENDIX E FRAME RATE MEASUREMENT CODE 


/*time.c****times digitize,calclsf and background, cetroid jitter*/ 
#define TITLE "Time calculations” 


#define AUTHOR "By W. J. Rall” 
/**~Quick C eet Oaipie of time COUEGRA ****** EARS See 
#include <time.h> 


while ( !kbhit() ) 

{ 
time1 = clock(); /* start clock */ 
for(i=O; i<128; i++) 


digitize sub frame(); 
calc_Isf(sub_frame); 
background (); 
centroid(); 


jit(i); 
time2=clock(); /* stop clock */ 


deltime = (time2-time1)/CLK_TCK; 
printf("\nTime for 128 frames =%f sec", deltime); 
} 
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